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Abstract. A well-known result with respect to the one dimensional nearest- 
neighbor symmetric simple exclusion process is the convergence to fractional 
Brownian motion with Hurst parameter 1/4, in the sense of finite-dimensional 
distributions, of the subdiffusively rescaled current across the origin, and the 
subdiffusively rescaled tagged particle position. 

The purpose of this note is to improve this convergence to a functional 
central limit theorem, with respect to the uniform topology, and so complete 
the solution to a conjecture in the literature with respect to simple exclusion 
processes. 



1. Introduction 

Informally, the one dimensional nearest-neighbor symmetric simple exclusion 
process follows a collection of random walks on the lattice Z which move inde- 
pendently except in that jumps to already occupied sites are suppressed. More 
precisely, the exclusion model is a Markov process r\ t = \r\t{x) : x S Z} evolving on 
the configuration space X = {0, 1} Z with generator, 



where rj x ' x+ is the configuration obtained from r\ by exchanging the values at x 
and x + 1 , 

!rj(z) when z ^ x,x + 1 
r)(x) when z — x + 1 
T](x + 1) when z — x. 
A more formal treatment can be found in Liggett [9]. Later, in Section 2, we will 
also give Harris's description of the model in terms of a "stirring process." 

As the process is "mass conservative," that is no birth or death, one expects 
a family of invariant measures corresponding to particle density. In fact, for each 
p G [0, 1], the product over Z of Bernoulli measures v p which independently puts 
a particle at locations x S Z with probability p, that is v p (j] x = 1) = 1 — v p (t\ x — 
0) = p, are invariant (cf. Liggett [5])- 

In this note, we concentrate on the integrated flux across the origin J(t), and the 
position of a tagged, or distinguished particle X(t), say the first particle to the left 
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of 1/2, when initially the exclusion process starts in an equilibrium v p for < p < 1. 
Both objects are interestingly connected, and have been long well-studied in the 
literature (see Section 8.4 in Liggett [S], Section 6.4 in Spohn [TJ], De Masi- Ferrari 

HI)- 

Perhaps the most intriguing behavior of the current and tagged particle is their 
subdiffusive fluctuation behavior, explained physically in part by the enforced or- 
dering of particles with no leapfrogging allowed in the dynamics. It was shown in 
Arratia [I], Rost-Vares [TJ] and De Masi- Ferrari [3] that 

t-^ 4 J{t) 4 N(Q,a 2 j) and t-^ 4 X{t) ^ N(0,a 2 x ) (1.1) 

where Oj = v/2/7r(l — p)p and a\ — y/2/n(l — p)p^ 1 - This is in contrast to 
the diffusive behavior in higher dimensions or when the jump probability is longer 
range (cf. Chapter 6 in Spohn [14], Part III in Liggett [10], Sethuraman [13]). 

Given these results, the general belief (see Conjecture 6.5 in Spohn [14]) is that 
the process limits with respect to the rescaled current and tagged particle position 
converge to respective fractional Brownian motions with Hurst parameter 1/4. In 
fact, straightforward modifications of the arguments for (|1.1[) give convergence in 
the sense of finite-dimensional distributions (a case of Theorem 1.2 Landim-Volchan 
[5] gives a specific statement; see also Jara-Landim [7]). However, it appears the 
full functional central limit theorem conjectured in Spohn [T3], with respect to 
exclusion processes, has not been addressed. 

The aim of this article is to complete the proof of this conjecture by supplying 
path tightness estimates to show, as A f oo, that both 

o-7 1 A- 1 / 4 J(Ai) => B 1/4 (i), and o-£\-V*X{M) => M lfi {t) (1.2) 

where IBi/4(£) is the standard fractional Brownian motion with parameter 1/4, and 
=> denotes weak convergence in -D([0, 1]) endowed with the uniform topology. 

The plan of the paper is to give some preliminary representations and estimates 
in Section [21 and then deduce the limits (|1.2|) . through certain maximal inequalities 
and discrete time process approximations, in Corollaries 13.51 and 14.31 in Sections [3] 
and [H respectively. 

Throughout, unless otherwise clear, P = P u and E — E v denote the process 
measure and expectation starting under the equilibrium v p . Also, as standard, \_x\ 
denotes the integer part of x. 

2. Representations of current and tagged particle 

In this section, we state a convenient construction of the exclusion process 
through a "stirring process" first introduced by Harris [6], and discuss some repre- 
sentations of the current and tagged particle position. 

2.1. Stirring process. The stirring process £J £ 1 for i e Z is defined as follows. 
At time t = 0, a particle is put at each site and we define £g = i for each i £ Z. 
To each bond (x, x + 1) with x E Z, we associate a Poisson process (clock) with 
parameter 1/2. When the clock rings at bond (x, x + 1), the particles at these sites 
interchange their positions. Then, Q is the position at time t of the particle which 
was at i at time 0. Given an initial configuration rj, the simple exclusion process, 
in terms of the stirring process, is 

Vt (x) = l{xe{%:r)(t) = l}}, 
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that is, in words, r]t(x) = 1 if and only if there is an i 6 Z so that £1 = x and 
17(0 = 1. 

2.2. Current representations. Let N + (t), N_(t) be counting processes, with in- 
finitessimal rates (1/2)?7 S (0)(1 — r) s (l)), (1/2)?7 S (1)(1 — ?y s (0)), which count the num- 
ber of particles which cross — » 1 and 1 — > up to time t respectively. Then, the 
current across the bond (0, 1) up to time t is given by 

J(t) = N + (t)-N^(t). 

As N+(t) - | X * 77 s (0)(l - r) a (l))ds and N-(t) - \ f* Tf a (l)(l - ?? s (0))ds arc martin- 
gales, we note the useful decomposition 

J(t) = M(t) + A(t) 

where M(t) = N+(t) - N-(t) - \ / *(r? s (0) - r] s (l))ds is a martingale, and A(t) = 

I Jo(%(°) - ^(!))^ (cf- (IH.2.37) in Liggett 

In terms of the stirring process, following the development in De Masi-Ferrari 
[4], define 

K + (t) = ^l{?j>0}; K-(t) = £l{g<0}. 

i<0 i>0 

In words, K + (t), K~(t) are the number of stirring particles starting to the left and 
right of the point 1/2 and sitting at the right and left of 1/2 at time t respectively. 
Denote 

n m - J ^ > °^ when * - 

' \ l{£t < 0} when i > 0. 
As in the stirring process all sites are always occupied, each crossing of the 
bond (0, 1) in one direction corresponds to a simultaneous crossing in the opposite 
direction. Then, K + (t) — K~(t) is constant in t, and since K + (Q) = K~(0) = 0, 
K+(t) = K~(t) := K(t), for all t > 0. Hence, 

m = E H£ > - E ^ 

For > 1, let ii < i 2 < ... < ix(t) < be the random locations for which 

> 0, and < ji < j% < ... < JkU) be- the random locations for which £;? fc < 0. 
Define = r)(ik) and B^ = rj(jk) and Ak{i) = Bt — B^ . Then, with the 
convention that the sum from 1 to is equal to 0, we have the representation 

K(t) 

j(t) = E Mt). 

fe=l 

We now state some known facts and consequences. 

(a) Clearly, given K(t), and the random locations {ik} and {jk} with 1 < k < 
K{t) the variables {Ak(t)} are independent, identically distributed, mean 0, and 

P(Ak(t) = 1) = P(A k (t) = -1) = p(l - p), and P(A k {t) = 0) = 1 - 2p(l - p). 

(b) Also K{t) is a sum of negatively correlated 0, 1 valued random variables. 
Moreover, by Lemma 4.12 Liggett [S], for all finite T C Z, and all AcZ, 

p(p|(£e^)) < U P (^ £A )- (2- 1 ) 
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(c) From basic considerations, E[K(t)] = _E(z(0,i)+), where z(0,t) is a symmet- 
ric random walk on Z starting at the origin. Then, E[K(t)] < \/t, and 

r E[K(t)] 1 
hm = — = — ==. 

t^o° y/t V27T 

(d) We have also some p-moment estimates. Denote ||V|| P = {E\V p ]) 1 / p for 
simplicity. 

Lemma 2.1. For integers p > 1, there is a constant Co = C'o(p) < oo so that for 
allt>l, 

\\K{t)\\ p < C Q Vt. 
Proof. First, as K(t) = J2i<o U%{t), and by (|2.1[) . we have that 



\K{t)\\ p < \ Y,Ui 



i<0 



r ~ T d 



where Ui = Ui(t) and {Ui} are 0, 1 valued independent random variables. 

To further estimate, we use Rosenthal inequality (Theorem 1.5.9 in De la Pena- 
Gine [3]) on independent and nonnegative random variables {Pi}: For p > 1, there 
exists a constant C\ — C\ (p) < oo such that 

Now, applying Rosenthal's inequality, for a positive integer p, noting that Uf = Ui, 
we have, 

llE^I ^ C 1 (j2m]+[E E P!}\ 1/P ) < C^ElK^ + iElKit)})^ 

i<0 P i<Q i<0 

The result follows from properties of E[K(t)) listed in part (c), and that p > 1. V 

2.3. Tagged particle representations. Consider a distinguished particle in the 
exclusion system. One representation for its displacement Z(t) = X(t) — X(0) is 
through the "Lagrangian frame," Ct = Ozfj^Vt where Q y r\ is the shifted configuration 
{6 y T]){x) — f] x +y Then, ( t is a Markov process on £' = {( e S : ((0) = 1} with 
generator 



(£0(0 = \ (HC' x+1 )-H0) 

xjt-1,0 

+ 1 E a- co (0(7*0-0(0) 



i=— 1,1 

where Tfc£ is the configuration obtained by displacing the tagged particle k steps 
and then shifting the frame, 

{((x + k) when i^O, — k 
C(0) when x = 
£(fc) when x = —k. 

Define Af + (t),Af-(t) as the counting processes, with infinitessimal rates (1/2)(1 — 
Cs(l))) (1/2) (1 — Cs( — which count the number of frame shifts of size 1 and —1 
respectively up to time t. Then, 

X(t)-X(0) = Af+(t) - Af-(t). 
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Similar to the current representation, Af+{t) — \ J* *(l — £ s (l))<is and J\f- (f) — | J (l — 
Cs( — l))ds are martingales, and 

X(t)-X(0) = M{t)+A{t) 

where M(t) = Af+(t) - W-(i) - § / * (&(-l) - C*(l))ds is a martingale and A{t) = 
\ J '(Cs(-l) ~ Cs(l))ds (cf. Proposition III.4.1 in Liggett [TD]). 

On the other hand, with respect to the stirring process and a configuration 77 
drawn from v p , following Diirr-Goldstein-Lebowitz [5] and the exposition in De 
Masi- Ferrari [4j , for k > 1 let Yfc (t) be the position of the feth particle of % to the 
right of 1/2; for k < let lfe(£) be the position of the (\k\ + l)th particle of rj t to the 
left of 1/2. Then, at time t, the tagged particle, initially the Oth labeled particle, 
is the J(i)th particle, 

X(t) = Y J(t) (t), 
where J(t) was defined in Subsection 2.2. 

It will also be useful to note, under the invariant measure i/ p , that Y n (t) = Y n (0), 

v m = / Yl{t) + di{t) for n - 2 

n[ ' I Y (t)-^72 n di(t) forn<-l 

where di(t) = Yi + i(t) — Yi is the spacing between the ith and (i + l)th particles, 
and also, {rfj(t) : i ^ 0}, and |5o(£)| + 1 have independent Geometric(p) 

distributions. 



3. Tightness and fBM limit for the current 

In this section, we prove the following theorem which is the main vehicle in the 
article. At the end of the section, we state as Corollary [33] the fractional Brownian 
motion invariance principle for the current. 

Theorem 3.1. Under initial distribution v p , the stochastic process \~ 1 / 4 J(\t) is 
tight in -D([0, 1]) endowed with the uniform topology. 

Proof. The proof follows from showing that the discretized version A -1 / 4 J(|_A£j ) is 
tight in D([0, 1]) in the uniform norm (Proposition ^. 3p . and then that the difference 
with the desired process is negligible (Proposition [33J ■ V 

The first step is to state a useful maximal inequality. 

Lemma 3.2. For integers m > 1, and even integers p > 2, 

E[J p (m)} < C 2 m p/i 

where C2 — C% (p) is a constant. Moreover, for even integers p > 6, there is a 
constant C3 = C^{p) such that 



E 



max J p (i) 

KKm 



< C 3 m p / 4 . 



Proof. We recall Marcinkiewicz inequality (Lemma 1.4.13 in De la Pena-Gine [3]): 
For p > 1, there exists a constant C4 = C±(jp) such that for centered, independent 
LP random variables {Pi}, 
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Denote by T the cr-algebra generated by K(t), and the random locations 
and {jk} with 1 < k < K(t) introduced in Subsection 2.2. Then, by conditioning 
first on taking into account that Ak(m) 2 < 1 for all k, m, and properties in part 
(a) Subsection 2.2, we have 

E[J p (m)} = E E ( Mm)] P 

l<k<K{m) 



T 

p/2- 



l<k<K(m) 

< C 4 E[K p ^(m)}. 

The first statement now follows by Lemma I2TT1 

For the maximal inequality, we note first, from stationarity of J(t) and the proven 
first inequality, for all integers < i < j < to, 

E[(J(3) - J(i)) p ] = E[j?(j - »)] < C 2 (j - i)*t\ (3.1) 
We now recall a case of Theorem 3.1 in Moricz-Scrfling-Stout [TT]: Let Sij = 
~Y^k=i$k where {/3k} are arbitrary random variables. Let also // > 1 and a > 1. 
Suppose for some nonnegative numbers {tife}, ElSij]^ < (J2l=i Uk ) a f° r au 1 — 
i < j < n. Then, there is a constant C5 = C^^/j,, a) such that 

, n x a 

E[xattK{\Si, 1 \,...,\Si, n \y t ] < C 5 (J2 Uk j ■ ( 3 - 2 ) 

^ k=l ' 

Then, as p/4 > 1, applying (|3.2| with respect to (|3.1|) . the second statement 
follows. V 

We now consider the discretized process. 

Proposition 3.3. Under initial distribution v p , the stochastic process A -1 / 4 J( |_AtJ ) 
is tight in -D([0, 1]) endowed with the uniform topology. 

Proof. According to Billingsley [2] , a well-known tightness condition is to show, for 
all e > 0, that 

limlimsupPf sup A~ 1/4 | J{[Xt\) ~ J{[Xs\)\ > e ) =0 
v s,*e[o,i] 7 
which reduces, in our situation of stationary increments, to proving 

limlimsupfT 1 / 3 ! sup A~ 1/4 | J{[Xs\)\ > e ) = 0. 

A^oo ^£[0,5] J 

By Chebychev's inequality, 

P\ sup A- 1 / 4 |J(LAsJ)| >e ) < e- e E\ sup A~ 6/4 | J 6 (LAsJ )| ) . 
\se[o,s] J \se[o,s] J 

Also, by the maximal inequality in Lemma 13.21 

El sup A- 6 / 4 |J 6 (LAsJ)| ) = E\ max A- 6 / 4 |J 6 ^)| ) < C 3 6 3/2 
\se[o,s] J \^i<i<LA5j J 

which is enough to conclude the proof. V 
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The difference between the discretized and desired process is handled as follows. 



Proposition 3.4. 

lim p( sup \- 1/4 \J(\t) - J([\t\)\ > s) = 0. 

A^oo \o<Kl / 

Proof. By stationary increments, noting (nt) — \nt\ < 1, 
P( sup \- 1/4 \J(\t) - J([\t\)\> e) 

^0<t<l ' 

< P\ sup A- 1/4 |J(t) - J(s)| > e ) 



0<s<t<A 
' t-S<l 



< 3Vp( sup \- 1/A \ J(t) - J(i)\ > e/3) 
i=0 V 4 <t< 4 +i 

= 3(LAJ+1)P( sup | J(t)\ > eA 1 / 4 /3 



0<t<l 



< 6£ 



sup I J(<)| 4 /( sup |J(i)| > e\ 1/4 /3) 

0<t<l ^0<*<l ' 



We now show that -E[sup 0<(<1 |J(i)| 4 ] < oo to deduce that the last quantity 
vanishes as A t oo. Indeed, from the decomposition J(t) — M(t) + A(t) where M{t) 
is a martingale and A(t) = (1/2) J '(^o(s) — Vi( s ))ds, as the integrand 1 770 — Vi\ 1S 
bounded by 1, we need only bound _E[sup 0<t<1 M 4 (t)]. By using Doob's inequality 
and simple computations, 



E 



sup M 4 (i)l < (4/3) 4 £[Af 4 (l)] < 4(4/3) 4 (£[J 4 (l)] + lY 
o<t<i - 1 v ' 



But, J(t) = N+(i) — N- (t) is the difference of two counting processes each with rates 
bounded by 1/2, and so by coupling with respect to dominating Poisson rate(l/2) 
processes, we have -E[J 4 (1)] < 00; this can also be seen from computing directly 

with J(i) = Ek=iMV- v 



Corollary 3.5. Under initial distribution v p , with respect to the uniform topology 
on -D([0, 1]), we have as A j 00. 

\- 1/4 J{\t) =► crjB 1/4 (i) 

where B 1 / 4 (t) is the fractional Brownian motion process with index 1/4. 

Proof. From Theorem l3.il we know that A _1 ' 4 J(Ai) is tight. Also, from the liter- 
ature (cf. Landim-Volchan [8]) the finite-dimensional distributions of any limit are 
Gaussian. Hence, if W(t) is a limit along a subsequence, this limit is a continuous 
Gaussian process. By Lemma 13.21 the sixth moment of A _1 / 4 J(Ai) is uniformly 
bounded. So, we have uniform integrability for first and second powers of the 
process (see also De Masi- Ferrari [4]), and therefore convergence of these moments. 
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To finish, we just have to compute the limit of covariances 
4cov(V 1 / 4 J(Xt), A~ 1/4 J(As)) 



= E A- 1 / 4 J(At) + A- 1/4 J(As) 
= E A-^J^At)] +E\\- 1 / 2 J 2 (Xt) 
+ 2cov(A- 1 / 4 J(Ai),A- 1/4 J(As)) 
Then, for t > s, recalling 



A" 1 / 4 J(Xt) - \- 1/4 J(Xs) 



-E 



lim cov(A- 1/4 J(At),A- 1 / 4 J(As)) = Ufi + - . 



4. Approximation and fBM limit for the tagged particle 

In this section, we approximate X~ X I 4 X{X£) by X' 1 ^ 4 p^ 1 J(Xt) (Proposition ^. 2|) 
by adapting part of the proof of Proposition 2.8 in Durr-Goldstein-Lebowitz [5]. 
Hence, as a process, X-^ 4 X(Xt) will converge to the same fractional Brownian 
motion limit as X^ 1 ^ 4 p^ 1 J(Xt) (Corollary [ 



Je,k(t) 



The first step is to approximate on the integers. For < e < 1, fc > 1 and t > 0, 
define 

e\J(t)\ for \J(t)\ >t x l s + k 
Zp-^t^z + k) for \J(t)\ <tV 8 + fr. 

Recall from Subsection 12.31 the tagged particle representation X(t) — Yj(t)- 
Proposition 4.1. For < e < 1, we have 



\Y J(t) (t)-p- l J{t)\ 

k^roo V tgf + Je,k(t) 



lim PI sup 



> 1 



0. 



Proof. First, we note, as in Durr-Goldstein-Lebowitz [5], on the event 

G e , m = {\Y n (t) - p^n\ < e\n\ for all |n| > m\ 

that — p n\ < 3p m when \n\ < m because particles cannot cross and so 

\Y n (t)\ < max{|F_ ro (i)|,|Y m (t)|} < m^ 1 + e) < 2p- 1 m. (4.1) 

By the representation from Subsection 2.3, and that {di(t) : i ^ 0}, Yi(t) and 
|y (i)| + 1 a re independent Geometric(p) random variables, we deduce 

P(\Y n {t) — P~ lr n\ > \n\e for some \n\ > m 

i-i 

< 2Y / P(\(Yi(t)-p- 1 )+Y / ( d *( t )- 

l>'rn i—1 

< C 6 m- q 

where Cq = Ce(e, p, q) for any power of q > 0. 



> le 



FBM LIMITS IN SSEP 



9 



Hence, noting the remark made near (14. 

p(for some t G Z+, |K J(t) (i) - p" 1 J(t)| > |J e , k (i)|) 
<2 ^ P(\Y n [t) - /9 _1 n| > e|n| for some |n| > i 1/s + 

With g > 8, the last expression vanishes as k t oo. V 
Proposition 4.2. For <5 > 0, we Ziawe 

lim p( sup A- 1 / 4 |X(At) -p" 1 ^)! ><s) =0. 

A-+oo Vfe[0,l] ' 

Proof. First, we write 

|AT(A*) - p-^CAt)! < |X(Ai)-X(LAtj)| + |X(LAtJ)-p- 1 J(LAtj)| 

+ p- 1 |J(LAtJ)- J(At)| 

and note both 



lim P 

A— »oo 



( sup \- 1/4 \X(\t) - X([\t\)\> 5/3) = 

V t6[0,l] 7 



and 

lim p( sup A _1 / 4 |J(**) - J(|_AiJ)| >8/3) = 0. 

A-+oo V te [ .l] / 

Indeed, the second limit is estimated in the proof of Proposition 13.41 and the first 
limit is similarly argued: Write X(t) - X(0) = M(t) + A(t) where Mit) is a 
martingale and A(t) is an additive functional with integrand bounded by 1/2; note 
|X(0)| + 1 is a Geometric(p) random variable, and so -E[JT 4 (0)] < Cj; then, 

E\ sup X 4 (t)] < 4E[ sup (X(t)~X(Q)) 4 ] + 4P[X 4 (0)] 

0<t<l 0<t<l 



< 16(4/3) 4 (P[M 4 (1)] +lj +4C 7 

< cJe[X 4 (1)]+i); 



and, as X(t) — X(0) = Af+{t) — A/l(i) is the difference of two counting processes 
whose infinitessimal rates are bounded by 1/2, by coupling with respect to domi- 
nating Poisson rate(l/2) processes, P[X 4 (l)j < 4E[(X(1) - AT(0)) 4 ] + 4C 7 < oo. 
Hence, we need only show 

lim p( sup A- 1 / 4 |X(LAtJ) -p _1 J(|At|)| > 5/3) = 0. 
a^oo H e [o,i] ' 

This limit is the same as 

£ := lim Pf max \- x ^\X{l) - p' 1 J(l)\ > 6/3). 

A^oo \0</<LAJ / 

As X(l) = Y J{l) (l), and 

A \Y J{l) (l)- P J(l)\ - -— 
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for < e < 1 and k > 1, noting for large A that Sp^A 1 / 8 + fc)A~ 1/4 < 5/6, we 
have 

£ < ....pf max i4#>^ + p(sup' Fj '''' t ;-;;' J( ' )l >iY 

" A >oo y 0<Z< LAJ A 1 / 4 3 J \ teZ+ Je,k(t) J 

With e fixed for the moment, £ is further bounded by 

r J ^ s\ , / |y J(t) (t)- P -V(t)| \ 

hm P sup > — + hm P sup — > 1 . 

A ^°° \te[o,i] A 1 '" 1 3ey fe^oo y t £Z + Je,k(t) J 

The second limit vanishes by Proposition 14. II Also, the first limit, by the invari- 
ance principle already proved for A -1 / 4 J(Xt) with respect to continuous fractional 
Brownian motion (Corollary [33]), vanishes by later taking e J, 0. V 



As mentioned in the beginning of the section, from Proposition ^. 21 and Corollary 
we obtain the fractional Brownian motion limit for the rescaled tagged motion. 



Corollary 4.3. Under initial distribution v p , 

\- x ^X{\t) a x M 1/4 (t) 

in -D([0, 1]) endowed with the uniform topology, where B 1 / 4 (i) is the fractional Brow- 
nian motion process with parameter 1/4. 
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